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We discuss theoretical and experimental preparations for an indirect new physics search using 
the rare decay Bj K*^\-^ K7i)jJ.^jJ.^ focusing on CP violating observables. The separation of 
new physics effects and hadronic uncertainties is the key issue when using flavour observables 
in a new-physics search. Our analysis is based on QCD factorization and soft-collinear effective 
theory and critically examines the new physics reach of those observables via a detailed error 
analysis due to scale dependences, form factors, and other input parameters; we also explore the 
experimental sensitivities at LHCb using a full-angular fit method; finally, we make the impact of 
the unknown A /mt, corrections manifest in our theoretical predictions. 
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1. Introduction 

At the beginning of the LHC era and close to the end of the B factories at SLAC [|l|] and 
at KEK ^ and of the Tevatron B physics experiments 0], all experimental data on flavour 
mixing and CP violating phenomena are consistent with the simple CKM theory of the standard 
model (SM) [^], which means that all flavour- violating processes between quarks are governed by 
a 3 X 3 unitarity matrix, usually referred to as Cabibbo-Kobayashi-Maskawa (CKM) matrix [^]. 
The CKM matrix is fully described by four real parameters, three rotation angles and one phase. 
It is this phase that represents the only source of CP violation in the SM and that now allows for 
an unified description of all the CP violating phenomena. This is an impressive success of the SM 
and the CKM theory and can be illustrated by the overconstrained triangles in the complex plane 
which reflect the unitarity of the CKM matrix. The successful theory was honored by last year's 
nobel prize in physics [^. Thus, the CKM mechanism is the dominating effect for CP violation 
and flavour mixing in the quark sector; however, there is still room for sizable new effects and new 
flavour structures because the flavour sector has only been tested at the 10% level especially in the 
b ^ s sector. In particular, CP violating observables are a good testing ground for new physics 
scenarios. While the SM is very predictive by describing all CP violating phenomena via one 
parameter, many new physics models offer many new CP phases. 

In Ref. |||, ^, we worked out the theoretical and experimental preparations for an indirect new 
physics (NP) search using the rare decay B^ — )• K*t^i^ based on the QCDf/SCET approach: QCD 
corrections are included at the next-to-leading order level and also the impact of the unknown A/rri}, 
corrections is made explicit. The full angular analysis of the decay B^ — )■ K*{-^ K^7i^)£^i^ at the 
LHCb experiment offers great opportunities for the new physics search. New observables can be 
designed to be sensitive to a specific kind of NP operator within the model-independent analysis 
using the effective field theory approach. The new observables Aj\ Aj \ and A^' are shown to 
be highly sensitive to right-handed currents. Moreover, it was shown that the previously discussed 
angular distribution A ^ ^ cannot be measured at either LHCb or at a Super-B factory. 

In the present letter we extend this preparation work to CP violating observables in the rare 
decay. We already anticipate here that, in contrast to claims in the literature, the new physics reach 
of such CP violating observables is rather limited. More details of our analysis with further results 
will be published in a forthcoming paper [p^. 



2. CP asymmetries 

The decay Bd K*^i+£- with Z*" K-71+ on the mass shell is completely described by 
four independent kinematic variables, the lepton-pair invariant mass squared, q^, and the three 
angles 61, 6k, <p- Summing over the spins of the final particles, the differential decay distribution 
of B4 -> K*^£+e^ can be written as 

TTi — ^ — 5— 7:i: = ^-^(^''^''^jf*''^)' ^2.1) 

aq^ a cos did cos Ok* a f 32n 
The dependence on the three angles can be made more explicit: 



2 



New physics reach of CP violating observables in — > K*t^t 



Tobias Hurth 



J{q^, di,dK,(l)) = sin^ + 7f cos2 Ok + (/^sin^ 0^ + 7|cos2 0^) cos 20/ +J3sm^ 0jf sin^ 0/ cos 20 
+74 sin 20a: sin 20/ cos 0+75 sin 20/f sin 0/ cos + (7^ sin^ Ok + 7^ cos^ Ok) cos 0/ 
+7? sin 20^ sin 0/ sin0 +78 sin 20/5- sin 20/ sin0 +79sin^ 0/s:sin^ 0/ sin 20 . (2.2) 

The angles are defined in the intervals 



1 ^ cos 0/ ^ 1 , 



-1 ^cos0^ ^ 1 



(2.3) 



where in particular it should be noted that the angle is signed. The corresponding decay rate for 
the CP conjugated decay mode B° — ;> K*'^{-^ /r+7r")/x+/i" is given by 



dq^ d cos did cos Ok* dij) 32k 



J{q\di,dK*,(l>). 



(2.4) 



As shown in |[TT[], the corresponding functions Ji{q^,di,dK, 0) are connected to functions 7/ in the 
following way: 

■/l,2,3,4,7 — ^ •/l,2,3,4,7, ■/5,6,8,9 " " -^S.e.S.g (2-5) 

where 7/ equals 7/ with all weak phases conjugated. 

The 7/ depend on products of the seven complex K* spin amplitudes, Aj^^/^, A\\l/r^ ^ql/r^ 
At with each of these a function of q^. The amplitudes are just linear combinations of the well- 
known helicity amplitudes describing the B Kn transition. They can be parameterised in terms 
of the seven B ^ K* form factors by means of a narrow-width approximation. They also depend 
on the short-distance Wilson coefficients C/ corresponding to the various operators of the effective 
electroweak Hamiltonian. The precise definitions of the form factors and of the effective operators 
are given in Refs. [|8|, 10]. Assuming only the three most important SM operators for this decay 



mode, namely O7, O9, and Oio, and the chirally flipped ones, being numerically relevant we have 

,2^ 



Al_L,R 
A\\L,R 



{{cf+cf')^{c\if+c\if')] 



-NV2{ml-ml.] 



10 



Mb + rriK* 



-cm} 



Ai(<?2) 



+ ^icf'-Cf')T2{q'] 
q^ 



Aql,r = -N / {IniK- \/ q^ 



{{cf-cf')^{cif-cf'')}>< 



X { (m| - m\* - q^ ) {niB + ntK* )A 1 (^^ ) - A A2 (^^ ) / {mb + m/f. ) } 



+2mt{cf' - Cf'){{mi + 3mj,. - q')T2{q' 



2 2 
nig — mj^. 



-Tsiq')} 



At = Nl"^/Vq^{2{C%f -C\{/')}A,{q^) 
where the Q denote the corresponding Wilson coeficients and 

A A A 00 00 

X =mB + mf^^.+q -2{mgmf^, +mf^*q +mgq ), 



(2.6) 



(2.7) 
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%,v*\W^^n\h-^-4- (2-8) 



\ 3-2i07r5,„3l 



In Refs. [11, 12], it was shown that eight CP- violating observables can be constructed by com- 
bining the differential decay rates of dT{B — > K^K^i^i^) and dT{B — > K^n^l^i^). Besides the 
CP asymmetry in the dilepton mass distribution, there are several CP violating observables in the 
angular distribution. The latter are sensitive to CP violating effects as differences between the 



angular coefficient functions, 7, —7,. As was discussed in Refs. [ |11| , |12[ ], and more recently in 



Ref. [13], those CP asymmetries are all very small in the SM; they originate from the small CP vi- 



olating imaginary part of A„ = (K/:)V,*.)/(VffoV,*.). This weak phase present in the Wilson coefficient 
Cg^^ is doubly-Cabbibo suppressed and further suppressed by the ratio of the Wilson coeficients 
(3Ci+C2)/C9f« 0.085. 

Another remark is that the CP assymmetries related to 75,6,8,9 can be extracted from {dV + 



dr) due to the property ([2.5|)., and thus can be determined for an untagged equal mixture of B 
and B mesons. This is important for the decay modes — ^ K*^{-^ K^nP)l^i'' and — ^ 0(— > 
K^K^)£^i^ but it is less relevant for the self-tagging mode Bd K*'-\-^ K+K^)i+i^. 

3. QCDf/SCET framework 

The up-to-date predictions of exclusive modes are based on QCD factorization (QCDf) and its 



quantum field theoretical formulation, soft-coUinear effective theory (SCET) []14|,|15|]. The crucial 
theoretical observation is that in the limit where the initial hadron is heavy and the final meson has 
a large energy [|l^] the hadronic form factors can be expanded in the small ratios Aqcd / "^fo and 
Aqcu/E, where E is the energy of the light meson. Neglecting corrections of order l/mt, and a^, 
the seven a priori independent B ^ K* form factors reduce to two universal form factors ^± and 
[ p^ p^ ]. These relations can be strictly derived within the QCDf and SCET approach and lead to 
simple factorization formulae for the B ^ K* form factors 

Fi = Hi^^ + (l)B(^Ti® (p^, + 0{A/mh) . (3. 1) 

There is also a similar factorization formula for the decay amplitudes. The rationale of such for- 
mulae is that perturbative hard kernels like //, or 7] can be separated from process-independent 
nonperturbative functions like form factors or light-cone wave functions 0,. 

However, in general we have no means to calculate A /mt corrections to the QCDf amplitudes 
so they are treated as unknown corrections. This leads to a large uncertainty of theoretical pre- 
dictions based on the QCDf/SCET which we will try to make manifest in our phenomenological 
analysis. 

The theoretical simplifications are restricted to the kinematic region in which the energy of 
the K* is of the order of the heavy quark mass, i.e. <^ m|. Moreover, the influences of very 
light resonances below 1 GeV^ question the QCD factorization results in that region. In addition, 
the longitudinal amplitude in the QCDF/SCET approach generates a logarithmic divergence in the 



limit ^ — )• indicating problems in the theoretical description below IGeV [14]. Thus, we will 



confine our analysis of all observables to the dilepton mass in the range, IGeV^ ^6GeV2. 
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Using the discussed simplifications the K* spin amplitudes at leading order in 1 /m^ and 
have a very simple form: 



VlNmeil-s] 

-VlNniBil - 
Nnig 



{{Cf + Cf] 



Irhh 



•eff 



c: 



-C;o)} + ^(cf 



At 



T(Cio- 

^{CyQ-c[Q)] + 2m{cf' 



c 



2/2 

q /m%, m 



■sY 



10 



^\\{Ek' 



(3.2) 



with s = q^/nig, rhj = mi/niB- Here we neglected terms of 0(m|.*). 

Most recently, a QCDf/SCET analysis of the angular CP violating observables, based on the 



NLO results in Ref. []14[], was presented for the first time []13[]. The NLO corrections are shown to 
be sizable. The crucial impact of the NLO analysis is that the scale dependence gets reduced to the 
10% level for most of the CP asymmetries. However, for some of them, which essentially start with 
a nontrivial NLO contribution, there is a significantly larger scale dependence. The ^^-integrated 
SM predictions are all shown to be below the 10^^ level due to the small weak phase as mentioned 
above. The uncertainties due to the form factors, the scale dependence, and the uncertainty due to 



CKM parameters are identified as the main sources of SM errors [ 13]. 



4. New physics reach 

The new physics sensitivity of CP violating observables in the mode — >• K*£'^£^ was dis- 



cussed in a model-independent way [ |13| ] and also in various popular concrete NP models |gOp. It 
was found that the NP contributions to the phases of the Wilson coefficients Cy, Cg, and Cio and of 
their chiral counterparts drastically enhance such CP violating observables, while presently most of 
those phases are very weakly constrained. It was claimed that these observables offer clean signals 
of NP contributions. 

However, we show that the NP reach of such observables can only be judged with a complete 
analysis of the theoretical and experimental uncertainties. All details of our analysis with further 



results will be published in a forthcoming paper []10|]. Here we will restrict ourselves to the most 



important issues. To the very detailed analyses in Refs. [ 13, 20] we add the following crucial points 



• We redefine the various CP asymmetries following the general method presented in our pre- 
vious paper [||] : An appropriate normalisation of the CP asymmetries almost eliminates any 
uncertainties due to the soft form factors which is one of the major sources of errors in the 
SM prediction. 

• We explore the effect of the possible A / mi, corrections and make the uncertainty due to those 
unknown A/mt, corrections manifest in our analysis within the SM and NP scenarios. 

• We investigate the experimental sensitivity of the angular CP assymmetries using a toy Monte 
Carlo model and estimate the statistical uncertainty of the observables with statistics correp- 
sonding to five years of nominal running at LHCb {lOfb^^) using a full angular fit method. 
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We discuss these issues by example the two angular asymmetries corresponding to the angular 
coefficient functions J(^s and J%; 

Afo = (/6s -/6.v) /d{T + t)/dq^, As = {h-h) /d{T + t)/dq^ (4.1) 

'Within the SM the first CP asymmetry related to /gs turns out to be the well-known forward- 



backward CP asymmetry which was proposed in Refs. [[1^, 19]. 

As a first step we redefine the two CP observables. We make sure that the form factor depen- 
dence cancels out at the LO level by using an appropriate normalisation: 



The 7, are bilinear- in the K* spin amplitudes, so it is clear from the LO formulae (2.6) that - 
following the strategy of Ref.|j8|] - any form factor dependence at this order cancels out in both 
observables. We note that has the same form factor dependence as 76.s but has larger absolute 
values over the dilepton mass spectrum that stabilizes the quantity. In Fig. |l| the uncertainty due 
to the form factor dependence is estimated in a conservative way (for more details see Ref. 



for A(,s defined in ( |4.1[ ) and for A^^ defined in ([4.2|). Comparing the plots, one sees that with the 



appropriate normalization this main source of hadronic uncertainties gets almost eliminated. The 
left-over uncertainty enters through the form factor dependence of the NLO contribution. Fig. ^ 
shows the analogous results for the observable Ag . 

In the second step we try to make the possible A /mi, corrections manifest in our final results. 
To explore the corresponding uncertainties we introduce a set of extra parameters for each spin 
ampUtude: 

A = Ai X (1 + Cie'^' ) + e'^A2 x ( 1 + Cjf'* ) , (4.3) 

where Ai 2 are the relevant sub-amplitudes and d is the weak phase. We assume that the subampli- 
tudes can each receive n A/mi, correction as well as an additional, currently unconstrained strong 
phase, 01 2. For the absolute size of the corrections we use a dimensional estimate fixing Ci,2 to 
be of order 5 — 10%. To access the effects of these uncertainties on the individual observables, we 
form an ensemble of theory predictions, where each amplitude is randomly assigned values of 0i 2 
and Ci,2 from a Uniform distribution over the specified ranges. It is assumed that the values of 
these parameters are not functions of q^. This ensemble is used to calculate a 66% confidence band 
for each observable by looking at the spread of predictions for each observable at each point in q^. 
The bands produced show the expected uncertainty on each observable given the estimated ranges 
for the unknown parameters. 

Within the SM, we have only one weak phase and the two subamplitudes are contructed in the 
following way; 

A = Asm(A„ = 0) X (1 +Cye'^') + (Asm(A„ / 0) -Asm(A„ = 0)) x (1 +C2e'^) (4.4) 

It turns out that in spite of this very conservative ansatz for the possible power corrections - we ne- 
glect for example any kind of correlations between such corrections in the various spin amplitudes 
- the impact of those corrections is smaller than the SM uncertainty in case of the two observ- 
ables Ag^ and Ag . In the left plot of Fig. |3| the SM error is given including uncertainties due to the 
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scale dependence and input parameters and the spurious error due to the form factors. In the right 
plot the estimated power corrections are given, which in case of the CP violating observable Ag^ 
are significantly smaller than the combined uncertainty due to scale and input parameters. Fig. ^ 
shows the same feature for the CP violating observable . This result is in contrast to the one 
for CP-averaged angular observables discussed in Ref. [|8|], where the estimated power corrections 
always represent the dominant error. As the reason for this specific feature one identifies the small- 
ness of the weak phase in the SM. Thus, one expects that the impact of power corrections will be 
significantly larger when NP scenarios with new CP phases are considered (see below). 

In the third step we consider various NP scenarios. Here we follow the model-independent 



constraints derived in Ref. [13] assuming only one NP Wilson coefficient being nonzero. We 



consider three different NP benchmarks of this kind: 

1. |CNP| =2. and0NP = ;j/8 

2. |CNP| = 1.5. and0NP = ;r/8 

3. ic'iol = 3. and ©'k, = 7r/8 

The absolute values are chosen in such a way that the model-independent analysis, assuming one 
nontrivial NP Wilson coefficient at a time, does not give any bound on the corresponding NP phase. 
Scenarios with larger phase values will be discussed in Ref. [Il0|]. Fig. § shows that the CP violating 
observable A^^ might separate a NP scenario (2), while the central values of scenarios (1) and (3) 
are very close to the SM. Moreover observable A^ seems to be suited to separate scenarios (1) and 
(3) from the SM. 

However, to judge the NP reach we need a complete error analysis within the three NP scenar- 
ios. Thus, let us consider the possible impact of unknown power corrections in these cases: With 
one new CP phase involved we work now with three weak subamplitudes in which possible power 
corrections are varied independently. 

The left plots in Figs. |6| and ^ show that the possible A /nib corrections have a much larger 
impact than in the SM and become the dominating theoretical uncertainty. However, the two CP 
violating observables could discriminate the specific NP scenarios with new CP phase of order n/% 
from the SM in view of the theoretical errors only. 

In the last step, we analyse the experimental sensitivity of the angular CP asymmetries using 
a toy Monte Carlo model. The right plots in Figs. ^ and ^ show the estimates of the statistical 
uncertainty of Ag^. and A\ with statistics corresponding to five years of nominal running at LHCb 
(10/^^^). The inner and outer bands correspond to la and 2a statistical errors. The plots show 
that all the NP benchmarks are within the 1 a range of the expected experimental error in case of 
the observable Ag^,, and within the 2a range of the experimental error in case of the observable Ag . 

Thus, our final conclusion is that while the prospects of NP discovery of the CP conserv- 
ing observable presented in ref. [j8|] both from the theoretical and experimental point of view are 
excellent, the possibility to disentangle different NP scenarios for the CP violating observables 
remains rather difficult. For the studied observables, LHCb has no real sensitivity for NP phases 
up to values of 7r/8 in the Wilson coefficients Cg, Cio, and their chiral counterparts via the rare 
decays B — )• K*t^i^ . Even Super-LHCb with 100/Z?^^ integrated luminosity does not improve the 
situation significantly. 
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Figure 4: SM uncertainty in A\ (left) and estimate of uncertainty due to A /mi, corrections (right, light grey 
(green) corresponds to C\^2 = 5%, dark grey (green) to Ci_2 = 10%). 




Figure 5: New physics scenarios, assuming one nontrivial NP Wilson coefficient at a time, next to SM 
prediction for Ag^^* (left) and Ag (right), for concrete values see text. 




1234561 23456 



Figure 6: A^f'^'- Estimate of uncertainty due to A /mi, corrections (left) and experimental uncertainty (right). 
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Figure 7: Ag : Estimate of uncertainty due to A /mi, corrections (left) and experimental uncertainty (right). 
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